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Abstract 

We have developed a new technique for photothermal displacement spectroscopy that is poten- 
tially orders of magnitude more sensitive than conventional methods. We use a single Fabry-Perot 
resonator to enhance both the intensity of the pump beam and the sensitivity of the probe beam. 
The result is an enhancement of the response of the instrument by a factor proportional to the 
square of the finesse of the cavity over conventional interferometric measurements. 

In this paper we present a description of the technique, and we discuss how the properties of 
thin films can be deduced from the photothermal response. As an example of the technique, we 
report a measurement of the thermal properties of a multilayer dielectric mirror similar to those 
used in interferometric gravitational wave detectors. 

PACS numbers: 
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I. INTRODUCTION 



Recently there has been much interest in under standing noise mechanisms in dielectric 
coatings of mirrors used in advanced interferometric gravitational wave detectors 
0,y,UQ|. One noise mechanism that could potentially limit the sensitivity, and hence the 
astrophysical reach, of an advanced detector is coating thermoelastic-damping noise 
0, llpj ]. In order to accurately predict the level of this noise and to design coatings in which 
this noise source is minimized, we need to have accurate values for the thermal properties 
of the thin films that make up these dielectric coatings. Specifically, we need to know the 
thermal expansion coefficient a and the thermal conductivity k of the coating materials in 
thin-film form. 

It is well known that the thermal conductivities of materials in thin films can differ 
markedly from those of the same materials in bulk form 
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. . . -j , and there is 

evidence of a similar deviation in the thermal expansion coefficient jl6, 17. - Thus, there 
is a need to directly measure the thermal properties of candidate coatings in order to select 
ones that exhibit the lowest thermoelastic-damping noise. 

We have constructed an apparatus to characterize dielectric coatings on mirrors for the 
purpose of selecting the best coating for an advanced interferometric gravitational wave 
detector. We use interferometric photothermal displacement spectroscopy, a technique well 
suited for measuring the thermal conductivity and expansion coefficient of a thin film [19I 

0, 0, HHQ, 

but our method differs from conventional photothermal-displacement- 
spectroscopy techniques in one important point. We use a Fabry- Perot cavity to substantially 
enhance both the pump-beam heating power applied to the sample and the sensitivity of 
the probe beam. 

In this paper we both describe the apparatus and discuss a simple and intuitive way of 
looking at the physics of the photothermal response over a broad range of frequencies. 



II. THE INSTRUMENT 



Traditional interferometric photothermal displacement spectroscopy uses separate pump 
and probe beams, sometimes with different wavelengths, to locally heat a sample and mea- 
sure the resulting thermal expansion. One widely-used configuration uses the sample as the 
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end mirror in one arm of a Michelson interferometer (see for example De Rosa et al. 
have observed the differential photothermal response in a pair of Fabry-Perot cavities, in a 
configuration originally conceived as an optical readout for a resonant-bar gravitational- wave 
detector. Their experiment nicely demonstrated the photothermal effect in Fabry-Perot cav- 
ities, verifying the theory of Cerdonio et al. 1(| for the frequency-dependent photothermal 
response of a homogeneous material. 

There is currently a need to characterize the thin-film dielectric coatings that are ex- 
pected to be used in advanced interferometric gravitational wave detectors. We have the 
need to measure both the thermal expansion coefficient and thermal conductivity of such 
coatings, and so we have constructed an instrument based on interferometric photothermal 
displacement spectroscopy to perform such measurements on candidate coatings. Our ap- 
paratus uses a single Fabry-Perot cavity, as shown in Figure [H with the sample forming one 



of the mirrors in the cavity 



26 



27]. 



Figure Q shows a diagram of our experimental apparatus. Both the pump and probe 
beams are provided by the same laser, so that both can resonate simultaneously inside the 
cavity. The pump and probe beams are distinguished from each other by having orthogonal 
polarizations. The first half- wave plate after the laser adjusts the polarization, determining 
how much goes into the pump or probe paths via a polarizing beamsplitter. We phase 
modulate the probe beam using an electro-optic modulator (E.O.M.) and lock it 
extension the pump beam) to the cavity by use of the Pound-Drever-Hall method 
We amplitude-modulate the pump beam using an acousto-optic modulator (A.O.M.) and 
use lock-in detection to measure the photothermal response from the error signal. 

For a pump beam sinusoidally-modulated at a frequency / that is much higher than the 
unity-gain frequency of the servo, the error signal is given by 



an d by 

mm- 



-8RP probe ^-J Q (P)Ji(P) |^sin27r/t| 



where R is the (RF) response of the photodiode, P pro be is the power in the probe beam, L 
and T are the length and finesse of the Fabry-Perot cavity, A is the wavelength of the laser 
light used, (3 is the phase modulation depth (in radians), and <5Lsin(27r/'t) is the deviation 



3 



Laser 



A/2 



Polarizing J 

i 
i 



Beam- 
splitter 



PUMP 



Faraday 
Isolator 



A.O.M. 
E.O.M. 



PROBE 



^7 



High-Frequency 
Phase Modulation 



Servo Filter 



©-- 



V 





Photo- 
detector 



Low-Frequency 
Amplitude Modulation 



SAMPLE 



Fabry- Perot 
Cavity 



Lock-In Amplifier 




Input 



L _ _ _ J 

Reference 



FIG. 1: A diagram of our experimental apparatus. 

in L from the ideal resonance condition 2L = NX. Jo and J\ denote Bessel functions of 
order zero and one. 

Thus the signal we are measuring is proportional to [T8L/X) , as opposed to simply (5L/X) 
for the conventional Michelson interferometer. The finesse T of the cavity, and hence the 
enhancement to the sensitivity to the photothermal displacement SL, can, in principle, be 
made quite large, with values of T = 10 5 achievable with standard techniques. 

In addition to enhancing the sensitivity to a given photothermal displacement 8L, the 
Fabry-Perot cavity also enhances the displacement itself. If P pU mp is the power in the pump 
beam, incident on the cavity, then the power incident on the sample inside the cavity is 
approximately Pi nc ~ TP vumv . As we shall see below, the photothermal response 5L is 
proportional to the total power absorbed by the sample, which in turn is proportional to 
Pi nc . Thus the Fabry-Perot cavity enhances the signal measured at the lock-in amplifier by 
a factor of T 2 . 
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This enhancement of the signal over conventional methods opens up many new possi- 
bilities for photothermal measurements. For example, it should be possible to measure the 
thermal properties of very-low-absorption samples that would otherwise be inaccessible to 
photothermal techniques. 

The sensitivity of this method is, in practice, limited by two things. First, it is limited by 
a combination of noise in the cavity length and laser-frequency noise. Intrinsic cavity length 
noise is thermal in origin and is typically small enough to be negligible for photothermal 



experiments except at the mechanical resonant frequencies of the cavity itself 30|. Extrinsic 
noise in the cavity can be suppressed by suitable isolation methods, including suspension 
from seismic isolation platforms and enclosure in a vacuum apparatus. 

A solid-state Nd:YAG laser, such as the Lightwave 126 31], typically has a frequ ency 
noise of Su(f) ~ lOOHz/VHz x (100Hz/ /), where / is the measurement frequency [32]. 
For a 30cm cavity, this corresponds to an equivalent length noise of 



5L eq _ freq (f) = ^Su(f) 



1 x 10- 13 4^L x 



100Hz 
~1~ 



V30cm J \1.06Afim J 



'Hz 

The fundamental shot noise limit for measuring the length of the cavity is much lower 
and, well below the cavity pole, is 



yhc \f\ 

O-'-'eq—shot ~ 



^ <F \J Pprobe 



8.1 x 10~ 19 ^= x 
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Second, the sensitivity is limited by the degree to which the pump and probe beams can 
be made orthogonal, and therefore isolated from each other. For the photothermal signals 
we have observed, we have had no difficulty in reducing this pump-probe cross-coupling noise 
to much less than the amplitude of our signal, even at the highest modulation frequencies 
and lowest signal levels. 

III. THE PHOTOTHERMAL RESPONSE 

A rigorous modeling of the photothermal response requires a complete, multidimensional 
solution to the diffusion equation to solve for the temperature distribution inside the sample, 
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then a solution to the Navier-Stokes equations to solve for the resulting thermal expansion. 



This has been done in a variety of contexts 



33, 



34, 
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, but the form of the solutions 



is somewhat involved, and fitting them to experimental data to extract a sample's thermal 
properties, especially those of a sample that includes a thin, inhomogeneous surface layer, 
is computationally intensive. 

We may gain considerable insight into the relationship between the photothermal response 
and thermal properties of a sample by making a simple estimate of the frequency depen- 
dence and magnitude of the photothermal response. If we apply a sinusoidal heat source to 
the surface of a material, thermal waves will propagate into the material, with a resulting 
sinusoidal thermal expansion over the volume being heated. A very general property of 
thermal waves is that they decay away as they propagate, with a characteristic decay length 
that is equal to the wavelength of the thermal wave. In this sense, they behave very much 
like electromagnetic waves propagating into a normal conductor 38[. For a homogeneous 
material with thermal conductivity k, mass density p, and specific heat C, this penetration 
depth l t is given by 

Thus we may assume, for a rough approximation, that the material only gets heated to a 
depth of if The thermal expansion of this part of the material is given, again approximately, 
by 

5L = £ t aAT. 

Let's consider a sinusoidally modulated pump beam with a power given by 

p _ p , p sin(27r/t) 

r pump r Q ~r j m ^ ' 

where P m is the peak-to-peak modulated power in the pump beam. (For 100% modulation, 
Pq — Pm/2, and P m is the total amplitude of the pump beam.) If P abs = WP m is the (peak- 
to-peak) power absorbed at the surface of the sample (W being the absorption coefficient), 
then the total temperature change AT over a single cycle is AT = (P a b s /2f)/(pCV), where 
V is the volume of the material getting heated. 

For high frequencies, where £ t <C r and r is the radius of the laser spot doing the 
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heating, we can approximate the volume by 



V ~ ixr 2 



which gives an approximate photothermal response 



a 1 1 



or 



5L hi _ f (f) ~ P abs -4, (1) 
where the characteristic frequency for this homogeneous substrate is defined as 

/- = -^h- (2) 

The high-frequency photothermal response is proportional to 1//, and its magnitude gives 
us the thermal expansion coefficient a, provided we know P a b s , p, and C. (Or it can give 
P a bs if we know a, p, and C, etc.) 

For low frequencies, we approximate the volume by Q| 

V ~ - f ~7r£ 

2 V3 ' 

which gives a photothermal response of 

en 

5L low -f(f) ~ Pabs"- (3) 

Note that the low-frequency response is essentially independent of frequency, with a natural 
"turning point" between the low- and high-frequency regimes at the crossover frequency f s . 
This turning point does not depend on P abs or a, is determined by the length scale r , and 
can be used to get the thermal conductivity k, provided we know r and pC. 

A careful derivation gives the complete response, valid for a homogeneous material at all 
frequencies Q, Q as 



, P a bs\ 1 a(l + cO 

oL{f) 



i r , r +o ° , mV 2 / 2 

aw / dv- 



(4) 



2 / 7T K 

where a denotes Poisson's ratio. 

For a material with a coating, there is an additional length scale, the coating thickness t. 
At sufficiently high frequencies the thermal penetration depth £ t becomes less than t, and the 
thermal waves essentially only sample the coating. When this happens, the photothermal 
response will be dominated by the coating and we have, approximately 

5L f^ft(f) ~ Pabs — ^J, 
K c J 
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where the subscript c denotes a coating property. The transition frequency that defines this 
regime occurs where it — t, or 



and the coating characteristic frequency is 



fc - rr^- (6) 



At low frequencies, where ^ 3> t, the volume of the sample that gets heated and expands 
is predominantly substrate, with the coating contributing very little to the overall pho- 
tothermal response. In this regime we expect the net response of the coated sample to be 
essentially indistinguishable from that of an uncoated sample. 

If t C r there will be two high-frequency regimes, as shown in Figure EJ one at 
moderately-high frequencies in which t <C it *C tq, where the photothermal response has a 
1/ / frequency response and is dominated by the substrate, and another at higher frequencies 
where it <^Lt r$. In this highest-frequency regime, the frequency dependence is again 1//, 
but the photothermal response is now dominated by the coating. The crossover frequency 
between these two regimes occurs the frequency where i t = t, or / = f t . 

We may interpolate between these two high-frequency regimes to obtain an approximate, 
high-frequency photothermal response of a coated sample. 



SL(f) » 



a(l + a)\ jj_ f a c (l + a c )\ f c 



f V «c J ft + f 



(7) 



This is the formula we will use to fit our high-frequency data to extract the coating thermal 
properties. 



IV. RESULTS 

We measured the photothermal response of two samples. The first was single-crystal, 
synthetic, c-axis sapphire. The second was identical to the first except for the addition of 
a 4/xm multilayer dielectric coating of alternating layers of SiC>2 and TiC>2. This coating 



formed a high-reflectivity mirror for infrared light with a wavelength of 1.064/zm 
samples, including the dielectric coating, were provided by CVI Laser, Inc. 39] 
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FIG. 2: The expected photothermal response of a coated sample. When the coating thickness 
is much less than the laser spot radius, there are three clearly defined frequency regimes. The 
magnitude of the photothermal response in each regime and the transition frequencies between 
regimes together allow us to determine the thermal properties of the sample. 

Both samples had a thin (200nm) layer of gold deposited on their surfaces. The purpose 
of this gold layer was threefold. First, this matched the optical reflectivities of the two 
samples so that all other aspects of the measurements would be the same. Second, the 
gold boosted the absorption for these initial measurements to enhance the photothermal 
response signal. And third, we wanted to ensure that all of the absorption occurred in the 
gold layer at the surfaces of the samples, as opposed to somewhere in the bulk or deep within 
the coating. Using the tabulated values of the thermal expansion coefficient, density, and 
specific heat for bulk gold |4fJ, and the reported thermal conductivity for a 200nm gold 
film we estimate that because the thermal conductivity of the gold layer is so good, 
and its thickness so small, its photothermal response will be negligible below 250MHz. At 
those frequencies, the observed response will be dominated by the substrate and dielectric 
coating. 

Agreement between theory and experimental data is excellent at all frequencies, giving 
us high confidence in the measurement techniques. Moreover, the values for a(l + cr) and k 
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FIG. 3: Photothermal response of both samples. Upper dots are experimental data for Sample 2, 
with the dielectric coating. The lowest theory curve is that of Reference 10( , which accounts for the 



substrate's photothermal response at all frequencies. The other two curves show high-frequency 
asymptotic behavior for samples with (Equation [7J) and without (Equation 0J) dielectric coatings. 
Noise in the data at high frequencies (low amplitudes) is due to optical crosstalk between the 
pump and probe beams. These effects are systematic and appear at the same levels and with the 
same frequency dependence in both samples. Gaussian error bars are smaller than the points as 
illustrated. 



are consistent with our expectations from tabulated bulk values. (See Table 1.) 



Figure El shows the photothermal response of both samples, where the upper set of data 
points are from the sample with the dielectric coating. The thermal properties used in 
this fit are given in Table 1. The two parameters we varied for the fit were the thermal 
conductivity k, which determined the rolloff frequency, and the product a(l + a), which 
determined the overall amplitude. The bulk value for the product pC was assumed for 
both coating and substrate, and the absorbed power was measured from the visibility of the 
cavity. We measure the sum of the transmitted and reflected power and subtract it from 
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the incident power to find the total power lost in the cavity. We then assume that this loss 
is due to absorption, rather than scattering. This assumption is only justified because of 
the relatively high-absorption gold film present in our sample. The response of the coated 
sample was essentially identical to that of the uncoated sample at low frequencies, which 
agrees well with our expectations that the dominant photothermal response in that regime 
is due to the substrate. This agreement also gives us confidence that the absorption is the 
same for the two samples, which was one of the goals in applying the thin gold coatings. 

At high frequencies, the response of the coated sample differs substantially from that of 
the uncoated one, and it agrees well with our expectations based on the thermal penetration 
depth as outlined in Section HTT1 Using the tabulated values of p and C for bulk Si0 2 and 
TiC>2 we extract an effective thermal expansion coefficient and thermal conductivity 
of a c (l + a c ) = (8.6 ± 0.6) x lO^- 1 and n c = (1.08 ± 0.15)W/mK. This value of the 
thermal conductivity is less than would be expected from the bulk values of k for SiC>2 and 
Ti02, which are 1.18W/mK and lOAW/mK, respectively, which together would give an 
effective thermal conductivity of 1.83W/mK. This reduction of the thermal conductivity is 
consistent with other observations, which indicate that the effective thermal conductivity of 



a material in thin-film form is often lower that its bulk value 
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3. 



Material 


p (kg/m 3 ) 


C (J/kgK) 


a (K- 1 ) 


k (W/mK) 


Si0 2 
Ti0 2 
Gold film 


2.2 x 10 3 [2] 
4.23 x 10 3 [40] 
19.3 x 10 3 


670 [2] 
688 [40] 
126 


5.5 x 10~ 7 [2] 
5 x 10" 5 [2] 


1.4 [2] 
10.4 [11] 
25 [11] 








a(l + a) (K- 1 ) 




Sapphire 
Coating 


4.0 x 10 3 [2] 
3 x 10 3 


790 [2] 
680 


(5.53 ±0.06) x 10" 6 * 
(8.6 ±0.6) x 10~ 6 * 


64 ± 2* 
1.08 ±0.15* 



Table 1: Relevant material properties used in this work. Entries with an asterisk (*) 

denote values derived from our data. 



V. CONCLUSIONS 

We have developed a new method of interferometric photothermal displacement spec- 
troscopy that is well-suited for studying thin dielectric films of the type used in optical 
coatings. Our method uses a Fabry-Perot cavity to enhance both the sensitivity of the in- 
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terferometric measurement and the amplitude of the photothermal response (by enhancing 
the pump-beam power). Together, these can potentially provide orders of magnitude more 
sensitivity than conventional photothermal displacement methods. In this paper we demon- 
strate both the technique and the data analysis by measuring the thermal conductivity and 
thermal expansion coefficient of a high-reflectivity Si0 2 — Ti0 2 film on a sapphire substrate. 
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